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Abstract — This paper concerns with the effect of radiative 
heat transfer on unsteady hydromagnetic free convective 
couette flow of a viscous incompressible electrically conducting 
fluid in presence of variable suction. Employing perturbation 
technique, the governing equations of the flow field are solved 
and the expressions for the velocity, temperature, skin friction 
and the heat flux i.e. the rate of heat transfer in terms of 
Nusselts number N„ are obtained. The effects of the important 
flow parameters such as radiation parameter F, magnetic 
parameter M, slip flow parameters h h lt 2 ; suction parameters 
oti, a 2 , Prandtl number P r etc. on the velocity and temperature 
of the flow field are analyzed and discussed graphically with the 
help of figures and tables. 

Index Terms- Convective, Couette, Heat source, 
Hydromagnetic, Suction. 

I. INTRODUCTION 

The phenomenon of hydromagnetic couette flow with heat 
transfer has been a subject of interest of many researchers 
because of its possible applications in many branches of 
science and technology. Flows through porous media have 
several engineering and geophysical applications such as, in 
the field of agricultural engineering to study the underground 
water resources; in petroleum industry to study the 
movement of natural gas, oil and water through oil channels 
and reservoirs; in the field of chemical engineering for 
filtration and purification processes. A series of 
investigations have been made by the researchers where the 
medium is either bounded by horizontal or vertical surfaces. 
Several researchers have analyzed such similar types of flows 
under various physical situations. Soundalgekar and Wavre 
[1] discussed the unsteady free convection flow past an 
infinite vertical plate with constant suction and mass transfer. 
Bejan and Khair [2] analyzed the heat and mass transfer by 
natural convection in a porous medium. Takhar et al. [3] 
studied the radiation effects on MHD free convection flow of 
a gas past a semi-infinite vertical plate. Attia [4] presented the 
transient MHD flow and heat transfer between two parallel 
plates with temperature dependent viscosity. Chamkha and 
his team [5] discussed the radiation effects on free convection 
flow past a semi-infinite vertical plate with mass transfer. 
Das and his co-workers [6] analyzed the hydromagnetic flow 
and heat transfer between two stretched/squeezed horizontal 
porous plates. Singh [7] studied the MHD free convection 
and mass transfer flow with heat source and thermal 
diffusion. 
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Nagraju et al. [8] discussed the simultaneous radiative and 
convective heat transfer in a variable porosity medium. Singh 
and Sharma [9] analyzed the MHD three dimensional 
Couette flow with transpiration cooling. 

The influence of moving magnetic field on three dimensional 
Couette flow was discussed by Singh [10]. Das and his team 
[1 1] analyzed the unsteady free convection MHD flow of a 
second order fluid between two heated vertical plates through 
a porous medium with mass transfer and internal heat 
generation. Makinde [12] showed the free convection flow 
with thermal radiation and mass transfer past a moving 
vertical porous plate. Singh et al. [13] described the MHD 
free convection transient flow through a porous medium in a 
vertical channel. Ogulu and Prakash [14] studied the heat 
transfer to unsteady magneto-hydrodynamic flow past an 
infinite vertical moving plate with variable suction. Das and 
his team [15] reported the unsteady free convective MHD 
flow and heat transfer of a second order fluid between two 
heated vertical plates through a porous medium. Das et al. 
[ 1 6] examined the effect of heat source and variable suction 
on unsteady viscous stratified flow past a vertical porous flat 
moving plate in the slip flow regime. Recently, Das and his 
co-workers [17] studied the unsteady transient MHD free 
convective mass transfer flow past an infinite vertical porous 
plate embedded in a porous medium in presence of suction 
and heat sink. 

The present study estimates the effect of radiative heat 
transfer on unsteady hydromagnetic convective couette flow 
of a viscous incompressible electrically conducting fluid in 
presence of variable suction and heat source. The effects of 
the important flow parameters on the velocity and 
temperature of the flow field are analyzed and discussed 
graphically with the help of figures and tables. 


II. FORMULATION OF THE PROBLEM 


Consider a two dimensional unsteady free convective 
magnetohydrodynamic flow of a viscous incompressible 
electrically conducting fluid between two vertical parallel 
porous plates placed at a distance h apart in the slip flow 
regime in presence of variable suction and radiative heat 
source. Let a time dependent suction 
v'{t')= —Vq ( l + sAe^‘ ) (1) 

be applied at the plate y=0 and the same injection velocity be 
applied at the plate y= 1 . We choose x-axis along the plate and 
y-axis normal to it. Under the above conditions the equations 
governing the flow are: 

Momentum equation: 


8u' 


- -v' 0 (l + sAe )^C *g/3(T'- n)+ v 
at ay dy 


d 2 u' <jB { 


(2) 
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Energy equation: 
dT' 
dt 


' V V P C p dy' 2 pC p dy' 


( 3 ) 


valid for an incompressible fluid. In order to solve Equations 
(7)-(8), we assume 


The boundary conditions of the problem are: 
, TT , dn' ST' q , 

1 dy' dy' k - 


u (y, t) = u 0 00 + sui 00 e m +• 
0(y, t) = 6>o 00 + s#i 00 «■“'+. 


(10) 

( 11 ) 


u'-u 2 =l 2 ^,t'=t;„ 

dy 


at v' = h , 


(4) 


where T/ =- ^—^-L , L being the mean free path and fUi 

Mi 

the Maxwell’s reflection coefficient. 

The radiative heat flux q r is given by 

dq r 


' 7 = 4{t' , 


dy 

where I = J k 


de 
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du 

dt 


(f + eAe ~ mt )— = G,.0 + ^ - Mu , 


dy 


dy- 


dt v dy P r dy 


( 7 ) 

(8) 


Using Equations (10)-(11) in Equations (7)-(8), we get the 
following zeroth order and first order equations: 

Zeroth order: 

du 0 _r n . 92 u 
--^—- ( ~r r Vo +— ~ 

dy dy 


— MUn , 


d0 o 1 0 2 O o 


■^0 


(5) 


is the absorption 


dy P r dy 

First order: 

du 0 duj 


( 12 ) 

(13) 


- cm, -A ^ - = G r e, + 

dy dy dy 


d 2 u ] 

— 2 Mu i > 


(14) 

(15) 


coefficient at the wall, P /,y is Planck’s function and 7, is the 

frequency, u ' is the velocity, T' is the temperature, B 0 is 
the uniform transverse magnetic field, P is the volumetric 
coefficient of expansion for heat transfer, P* is the 
volumetric coefficient of expansion for mass transfer, k is the 
thermal conductivity, v is the kinematic viscosity, C p is the 
specific heat at constant pressure, a is the electrical 
conductivity, g is 

the acceleration due to gravity, A is a real positive constant, t 
is the time and s is a small positive number such that sA«l. 
Introducing the following non-dimensional variables and 
parameters, 


dy dy P r dy 

The corresponding boundary conditions are 

. du 0 duj d9 0 r d0, 

u 0 =a 1 +h 1 ,u I =h I , = — 7, =0 atv = 0 

dv dy dy dy 

du n da, 

u n =a 2 +h 7 ,u 1 =h 2 , 0 o =O,0 I =O at y = l (16) 

dy dy 

The solutions of Equations (12)-(15) under boundary 
condition (16) are given by 

u{y,t)=[A 5 e m5y +A 6 e m6y -A u e miy -A l2 e m2y )+se- at (A 7 e m7y +A s e’" 8y 
+B,e miy +B 2 e m2y -B 3 e m}y -B 4 e m4y -B s e msy -B 6 e"« y ) (17) 

0{y,t)=( +A 2 d n2y )+se- at (A 3 e m3y +A 4 e'" 4y - A 9 e m/y -A 1( f 2y ) 

(18) 

The wall shear stress i.e. the skin friction at the wall is given 


T = 


^ du^ 


9 y) v 


f du, ^ 


y = 0 V y J y=0 


+ ee 


f du ^ 


v v y= 0 


(19) 


in Equations (2)-(3), we get the following non-dimensional 
equations 


Using Equation (17) in Equation (19), it is given by 
r =(m 3 A 5 +m 6 A 6 —m 1 A I1 -m 2 A 12 )+se^ al (m 7 A 7 +m s A 8 

+ m 1 B 1 + m 2 B 2 - m 3 B 3 - m 4 B 4 - m 3 B 3 -m 6 B 6 ) . (20) 

The rate of heat transfer i.e. the heat flux at the wall is given 
by 


N„ = 


00 

v 5 tL 0 


+ ee 

V ) v =n l 


( 21 ) 


where M is the magnetic parameter, G, is the Grashof number 
for heat transfer, F is the radiation parameter, P, is the Prandtl 
number, a.i and a 2 are the suction parameters and h\ and h 2 
are the slip flow parameters. 

The corresponding boundary conditions are: 
dudO 

u = a,+n , — , = -7 at v = 0 , 

1 1 dy dy 

u = a , H , 0 = 0 at y = h. (9) 

dy 

III. METHOD OF SOLUTION 

We now seek the solutions for Equations (7)-(8) under 
boundary condition (9) for a particular case R=l, which is 


'}’=u V y y=0 V s y y=0 

Using Equation (18) in Equation (21), it is given by 

N u = ( nijAj + m 2 A 2 ) + ee~‘ a (m 3 A 3 + m 4 A 4 - mjA 9 - m 2 A I0 ), (22) 

where 

P„ 

_i_ 


i,=~^f + ^P, 2 +4P r F 


= — '---Jp 2 +4 P r F , 
2 2 y 


^f + -JP r 2 - 4P r (co- F] , m 4 = - 4P r [co- F ) , 


; =-- + -yll + 4M , m 6 =----yll+4M , 

2 2 2 2 

= - 1 2 + 1 2 4 T +W Z • m s =~-^l + 4(M-a>) , 


7 7 


? m 2 
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(in-,e mi - nije " 12 ) 


A 2 =- 


Aje " ll 
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„ _ G, A, 


B, = 
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,C 1 = l — h I m 5 , C 2 =l — h 1 m 6 , 


+ m 6 - (M - co) 

C 3 =a 1 +A II (l-h I m I )+A I2 {l — h 1 m 2 ) , 

C 4 = e"' 5 ( 1 - h 2 m 5 ) , C 5 = e'” 6 (l - h 2 m 6 ) 

C 6 = a 2 + A n e m} {l - h 2 m I )+ A 12 e m4 {l - h 2 m 4 ) 

C y — 1 h i^yiy , C g — 1 hjlYlg , 

C 9 = -Bj(l - hjtrij)- B 2 (l ~ h 1 m 2 )+ B 3 (l -hjm 3 ) 

+ B 4 (l- h,m 4 )+B s (l- h 1 m 5 )+B 6 (l - h 1 m 6 ) , 

C 10 = e m? (l - h 2 m 7 ), C n = e ms [l - h 2 m 8 ), 

C 12 =—B I e mi (l—h 2 m 1 )—B 2 e m2 (l—h 2 m 2 )+B 5 e" ,3 (l—h 2 m 3 ) 

+ B 4 e m4 ( 1 -h 2 m 4 )+ B s e"' 5 ( 1 - h 2 m 5 )+ B 6 e m<5 ( 1 - h 2 m 6 ) . 

(23) 


IV. RESULTS AND DISCUSSIONS 


points. The magnetic parameter affects the velocity of the 
flow field to an appreciable extent. Figure 2 depicts the effect 
of magnetic parameter M on the velocity of the flow field. 
The magnetic parameter is reported to decelerate the velocity 
of the flow field at all points due to the magnetic pull of the 
Lorentz force acting on the flow field. Figures 3 and 4 
estimate the effect of slip flow parameters on the velocity of 
the flow field. Comparing the curves of both the figures, it is 
observed that the effect of growing slip flow parameters h\ 
and h 2 is to enhance the velocity of the flow field at all points. 
Figures 5 and 6 point out the effects of suction parameters a t 
and a 2 on the velocity field. Both the parameters u t and a 2 
lead to accelerate the velocity of the flow field at all points. 



Figure 1. Velocity profiles against y for different values of F 
with M= 1, G>= 1, /? y=0. 1 , h 2 =0A, c/.r(). 1 , a 2 = 0.2, P,A)J \ , 
A= 0.5, ©=0.1,t=0.1 and £=0.01 


The present study considers the hydromagnetic convective 
couette flow with heat transfer in the slip flow regime in 
presence of variable suction and radiative heat source. The 
governing equations of the flow field are solved for velocity, 
temperature, skin friction and heat flux and the effects of the 
flow parameter such as radiation parameter F, magnetic 
parameter M, slip flow parameters h \ and h 2 , suction 
parameters a t and a 2 and the Prandtl number P, on the 
velocity, temperature and skin friction are analyzed with the 
aid of velocity profiles shown in Figures 1-6, temperature 
profiles shown in Figures 7-8 and Tables 1-2 respectively. 

4.1. Velocity field 

There is a drastic change in the magnitude of velocity of the 
flow field with the variation of radiation parameter F, 
magnetic parameter M, slip flow parameters h u h 2 and 
suction parameters a. / and a? in the flow field. Figures 1-6 
clearly show the variations in the velocity field due to the 
change in the above parameters. 

Figure 1 elucidates the effect of radiation parameter F on the 
velocity of the flow field. An increase in radiation parameter 
is reported to enhance the velocity of the flow field at all 



y 


Figure 2. Velocity profiles againsty for different values ofM 
with G,.=l, F=().l , /?/=(). 1 , /i 2 =0.1, c/.r(). I , a 2 = 0.2, P r =0.71, 
A= 0.5, ©=0.1, t=0.1and s=0.01 
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y 


Figure 3. Velocity profiles against y for different values of /? , 
with M= 1, G r = 1, FMl. l , h 2 = 0.1, a t =0.1, a 2 = 0.2, P,.=0.71, 
A= 0.5, ®=0.1, 1=0. land s=0.01 



Figure 5. Velocity profiles against y for different values of a t 
with M= 1, G r = 1, F=0.1, h,= 0.1, /? 2 =0.1, a 2 =0.2, P,.=0.71, 
A=0.5, o=0. 1, 1=0. land e=0.01 



T 


Figure 4. Velocity profiles against y for different values of/? 2 
with M= 1, G,= 1, F=0.1, h ?=0. 1 , «/=0.1, «2=0.2, P~Q.1 \ , 
A=0.5, o=0.1, 1=0. land s=0.01 

4.2. Temperature field 

The major changes in the temperature profiles of the flow 
field are due to the variations in the radiation parameter and 
Prandtl number in the flow field. These changes in the 
temperature field are analyzed in Figures 6 and 7 
respectively. A close observation of the curves of both the 
figures shows that a growing Prandtl number or radiation 
parameter decreases the temperature of the flow field at all 
points. But with growing values of Prandtl number and the 
radiation parameter, the decrease in temperature is more 
significant in case of Prandtl number. 



Figure 6. Velocity profiles against y for different values of a 2 
with M= 1, G r = 1, F=0.1, /?,= 0.1, h 2 = 0.1, a,=0.1, P,~0.1l, 
A= 0.5, ©=0.1, 1=0. land s=0.01 

4.3. Skin friction 

Tables 1 and 2 elucidate the variations in the values of skin 
friction at the wall against magnetic parameter M with the 
variation of a j and a 2 respectively. A close observation of the 
values of Tables 1 and 2 depicts that a growing suction 
parameter ai decreases the skin friction at the wall for a given 
value of magnetic parameter M, while in case of the other 
suction parameter a 2 the effect reverses. It is further observed 
that a growing magnetic parameter decreases the skin friction 
at the wall at all points of the flow field for a given value of ai 
or a 2 . 
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Figure 7. Temperature profiles against y for different values 
of F with M= 1, G>= 1, h]= 0.1, h 2 = 0.1, ai=0A, a 2 = 0.2, 
P r =0.7l, .4=0.5, ®=0.1, 1=0. 1 and s=0.01 



Figure 8. Temperature profiles against y for different values 
of P r with M= 1 , G r =l,F=0.\,hi=0.l,h 2 =0.l, «/=0.1, a 2 = 0.2, 
A= 0.5, o=0.1, <=0.1 and £=0.01 

Table 1. Values of skin friction i at the wall against M for 
different values of «/ with G r = 1 , F=0. I , .4=0.5, a 2 =0.2, 
/?i=0.1, h 2 = 0.1, P r = 0.71, o=0.1, <=0.1 and £=0.01 


T 

M 

ai=0 

a i=0.1 

ai=0.3 

ai=0.5 

0 

0.31603 

0.04471 

-0.49790 

-1.04052 

0.5 

0.30734 

0.03208 

-0.51843 

-1.06895 

1 

0.29913 

0.02004 

-0.53814 

-1.09632 

3 

0.27030 

-0.02304 

-0.60974 

-1.19644 

5 

0.24667 

-0.05949 

-0.67181 

-1.28413 


Table 2. Values of skin friction x at the wall against M for 
different values of a 2 with G r = 1, F=0.1, A= 0.5, a/=0.1, 
h 7 = 0 . 1 , h 2 =0. 1 , T > ,.=0.7 1 , o=0.1, <=0.1 and £=0.01 


T 

M 

a 2 =0 

a 2 =0.1 

a 2 =0.3 

a 2 =0.5 

0 

-0.04887 

-0.00208 

0.09151 

0.18510 

0.5 

-0.05701 

-0.01246 

0.07663 

0.16572 

1 

-0.06484 

-0.0224 

0.06248 

0.14735 

3 

-0.09347 

-0.05825 

0.01217 

0.08259 

5 

-0.11852 

-0.08912 

-0.02996 

0.02906 


V. CONCLUSIONS 

We report the following results of physical interest on the 

velocity, temperature and skin friction at the wall of the flow 

field from the present study. 

1. An increase in magnetic parameter M decreases the 
velocity and also the skin friction at all points of the flow 
field. 

2. A growing radiation parameter F enhances the velocity of 
the flow field at all points. 

3. The effect of growing slip flow parameters h\ and h 2 is to 

enhance the velocity of the flow field at all points. 

4. The suction parameters a t and a 2 have an accelerating 
effect on the velocity of the flow field at all points. 

5. The effect of increasing Prandtl number P r or radiation 
parameter F is to decrease the temperature of the flow 
field at all points. 

6. A growing suction parameter ai decreases the skin 
friction at the wall for a given value of magnetic 
parameter M, while in case of the other suction parameter 
a 2 the effect reverses. Further an increase in magnetic 
parameter decreases the skin friction at the wall at all 
points of the flow field for a given value of ai or a 2 . 
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